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Abstract 

The spectral determinant D(E) of the quartic oscillator is known to satisfy a 
functional equation. This is mapped onto the ^-related F-system emerging 
in the treatment of a certain perturbed conformal held theory, allowing us 
to give an alternative integral expression for D(E). Generalising this result, 
we conjecture a relationship between the x 2M anharmonic oscillators and the 
A2M-1 TBA systems. Finally, spectral determinants for general \x\ a potentials 
are mapped onto the solutions of nonlinear integral equations associated with 
the (twisted) XXZ and sine-Gordon models. 



1 e-mail: p.e.dorey@durham.ac.uk 
2 e-mail: tateo@wins.uva.nl 



Since the discovery of quantum mechanics, the spectral problem associated with 
the homogeneous Schrodinger operator 



AM*) = (-^ + z 2M ) Mx) = E k Mx) (1) 

on the real line has been the subject of much attention, with a supply of papers 
which continues to this day: refs. [fjldlOfl offer just a small sample of this work. Given 
the apparent simplicity of the system, it is at first surprising that much of the most 
remarkable progress has been made relatively recently. In the following we will be 
guided by the theory developed by Andre Voros in , and we refer the reader to 
these articles for a detailed explanation of the subject. Here we summarize a few 
facts that will be needed later. The confining nature of the potential in (||) means 
that the spectrum {Ei} of the theory is discrete. The properties of this spectrum 
can be encoded into spectral functions, the simplest example being the spectral 
determinant 

00 / E \ 

D M (E)=D M (0)l[(l + — ) (2) 

The constant D M (0) = sin(7r/(2M+2))- 1 reflects the definition of D M zeta- 
regularized functional determinant (see @). Dm{E) is an entire function of E 
and the positions of its zeroes coincide, by definition, with the negated discrete 
eigenvalues of eq. ([[]). Despite the absence of any closed expression for the E^, 
precise information about the spectrum can be obtained by various means. The 
particular aspect that will be important for us is the fact that the functions Dm{E) 
satisfy certain functional equations similar to those previously obtained for 

related Stokes multipliers [ 1 1 1 . These were used in [R-H to derive sum rules relating 



the different eigenvalues, but their utility was limited by the difficulty in finding 
solutions to the equations within the class of entire functions. In this paper we 
point out a surprising link between these functional equations and other systems 
of equations which have arisen in the last few years in a very different context, 
namely the finite-size spectra of integrable 1+1 dimensional quantum field theories. 
Numerical work confirms the match, and we feel that this unexpected connection 
between two a priori disconnected topics deserves to be understood at a deeper 
level. 

We begin by reviewing some basic properties of the spectral determinants. From 
the Bohr-Sommerfeld approximation one can deduce the asymptotic positions of 
the zeroes E = —E^ : 

6 (^fc) M ~2vr(A; + l/2) , k -> 00 (3) 
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where /i = (M+1)/2M and 



1V1 1 l 2 ^ 2M / 



In addition, Dm(E) admits a semiclassical asymptotic expansion for \E\ — > 00 with 
I arg.E| < 7r-5, <5 > 0: 

00 , 
In ~ V a^ 1 - 2 ^ ; a = - . ° . (5) 

i=o VP 7 

Now suppose that M = 2. In this case D{E) = D<z(E) satisfies the following 
functional relation [m]: 

D(Ej- 1 )D(E)D(jE) = D(Ej- 1 ) + £>(£) + £>(j#) + 2 (6) 

where j = e 2l7r / 3 . Together with the asymptotics just described, this is strongly 
reminiscent of the properties of solutions to thermodynamic Bethe ansatz (TBA) 
equations |12|, 13||. Consider, for example, the perturbation of a theory of para- 



fermions by the thermal operator of conformal dimensions A = A = 2/(h+2). This 
results in an integrable massive quantum field theory, associated with the A^_i Lie 
algebra. There are h—1 particle species, with masses M a = Mi sin(jra/h)/ sin(n/h), 
a = 1 . . . h—1. Species a and h—a are charge-conjugate: a = h—a. The scattering 
theory is factorisable, with two particle S-matrix elements J14[ |: 



a+6-1 



S ab = [] W ' a,b = l...h-l (7) 



\a-b\ + l 
step 2 



where, in the notation of JT|, M = (p-l)(p+l), (p) = sinh(f-Mff )/sinh(f-i7^). 
Non-perturbative information concerning the finite-size scaling functions of the 
model in a cylinder geometry can be obtained using the thermodynamic Bethe 
ansatz technique [12|, I^,[l6|,|l7]]. The simplest instance [12,13|| expresses the ground- 



state energy E(M\,R) as — ttc(MiR)/6R , where 

h-l 



o " * poo 

c(r ) = — / ^ «W cosh 0L a (0) , (8) 

^ a=l ^ 



L a (0) = ln(l+e~ ea(0) ) , r = MtR, and m a = M a /M x . The functions e a (0), a 
1 . . . h—1 (known as pseudo-energies) solve the following equations: 



. ft.— 1 

£q (#) = m a r cosh0- —^</> afe *L a (0) (9) 



2 



with (t> ah (6) = -id 9 lnS ab (9) and g*f(0) = !°° 00 d9'g{9 - 9>)f(9'). Now con- 
sider Y a {9) = e £ °( fl ). These are entire functions of 9, with periodicity Y a {6 + 
in(h+2)/h) = Y a (9) |i~8| . Conjugation symmetry of the ground-state equations 
means that e a {9) = e a (9), so the Y's are entire functions of t = exp(2h0 / (h+2)) 
on the punctured f-plane C* = C \ {0}. In fact the Y's are thought (cf. jOJ) 
to be analytic functions of the variables a± = (r exp(±9)) 2h ^ h+2 \ with a finite 
domain of convergence about the point (a+,a_) = (0,0). The domain of con- 
vergence is finite because of square-root singularities linking the ground state to 
excited states [||,|I|,|l7|. It was also shown in |18| that the Y's satisfy a set of 
functional identities, known as a Y-system. At h = 4, t = exp(4#/3), and, taking 
the conjugation symmetry into account, the Y-system is: 

Yi(e- i7r/3 t)Y(e i7r/3 t) = 1 + Y 2 (t) (10) 
Y 2 (e- i ^H)Y 2 {e i ^h) = (l + Y^i)) 2 (11) 



Substituting ( ^0| ) into (pi]), it is easy to see that Y\(i) satisfies a constraint involving 
itself alone: 

Y 1 (e- 27ri / 3 t)Y 1 (t)Y 1 (e 27ri / 3 t) = Y^e-^H) + Y(t) + Yi(e 2 ^ /3 t) + 2 . (12) 

The relation with equation (||) is clear, but the analytic properties of Y\ don't quite 
match those of D yet. In particular, Y% has an essential singularity at t = 0. To 
remedy this, we take a massless limit, replacing the driving term m a r cosh 9 with 
m a re e (this amounts to setting a_ = 0). The resulting Y's are now nonsingular 
at t = 0, and furthermore, for the ground state their zeroes lie on the line Qro 9 = 
37r/4, the negative real axis in the t-plane. Setting rri2r = 6o|m=2 and t = E, 
and identifying Y\(t) with D{E), all of the standard properties of the spectral 
determinant of the quartic oscillator are reproduced. For example, the large 9 
asymptotic Y 2 {9) ~ b$e e is obtained by dropping the convolution term in @, and 
implies that Y 2 (9) takes the value —1 at the points 9 = + vri/2, where 

b e Xk ~ 2ir(k + 1/2) , k -> oo . (13) 

Combined with ([ToD, this shows that the zeroes of Yi(0) are at 9 = Xk + 37ri/4, 
matching the asymptotic behaviour (0). Finally, at t = the solutions of the Y- 
system are Yl = 2, Y 2 = 3, matching the result D(0) = 2. Still unsatisfied, we 
performed a numerical check. Eq. (H) was solved for real and then, as in [19|, 



eq. (g) and the Y-system were used to obtain the values of Y\(9) on the line 
^sm9 = 37r/4. The first zeroes were found to high accuracy, and the resulting 
predictions for the first five energy levels of the x 4 potential are compared with 
earlier results in Table I. 
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E k (TBA) E k (QM) 



1.06036209048418 1.06036209048418289965 a 

1 3.79967302980139 3.79967302980 6 

2 7.45569793798672 7.45569793798673839216° 

3 11.64474551137815 11.6447455114 6 

4 16.26182601885024 16.26182601885022593789 a 

5 21.23837291823595 21.2383729182 6 

Table I: Energy levels for the x 4 potential from the TBA, com- 
pared with previous results: a from |^|,^| , b from 

For M > 2 the equations satisfied by Dm(E) become more intricate, and we 
have yet to map them explicitly into known TBA systems. Instead we took a short- 
cut, though later we shall give an alternative, and more systematic, treatment of 
the problem. The functional relations for Dm{E) have a Z/j + 2 symmetry [|7j], where 
h=2M. This suggests an examination of Y-systems which share this symmetry 
in order to find a generalisation of the M=2 result. Of the diagonal scattering 
theories, this picks out the models associated with the A^—i or Lie algebras 

(see [1^, IE, 18 1), for which the Y-systems are 



Y a (9 - if)Y a {9 + il) = J] (1 + Y b {9)) lah (14) 



6=1 



where r is the rank and l a b the incidence matrix of the relevant Dynkin diagram. 
However, the constants Y a {6=— oo) do not match the value sin(7r/(2M+2)) _1 of 
Dm (0). But all is not lost: we can invoke another system of functional relations, 
related to the Y-system, called the T-system (cf. pC|): 

r 

T a (9 - q)T a (0 +q) = 1 + J] T h (9) la \ (15) 

6=1 

with Y a (9) = ilb=i T b(9) lab - When M=2 we have T 2 {9) = Y 1 (9), and so we can 
also search for our generalisation amongst the T-systems. Asymptotic checks lead 
to the conjecture that Dm{E) coincides with the function Tm{9) of the massless 
A2M-1 TBA system obtained from eq. (||) by setting h=2M, replacing the terms 
m a r cosh# by m a re e , and setting mMf = bo and = E. This was checked using 
the fact that the zeroes of Tm{9) on the line Qm9 = (h+2)ir/2h correspond to 
zeroes of 1+Ym(9) on the line Qm9 = tt/2, and these can be located using (||) and 
the Y-system as before. Tables II and III show some results for M = 3 and 4. 



4 



E k (TBA) 



E k (QM) 



1.144802453797075 

1 4.338598711513990 

2 9.073084560921449 

3 14.93516963491078 



1. 14480245379707° 
4.3385987115 6 
9.07309 c 
14.9351696349 6 



Table II: Energy levels for the x 6 potential from the TBA com- 
pared with previous results: a from table 2 of EO], b from table I 
of i.andc from table VII of B , rescaled by 2^/ 4 



k 


E k (TBA) 


E k (QM) 





1.2258201138005 


1.22582011382" 


1 


4.7558744139607 


4.7558 c 


2 


10.2449469772369 


10.2450 c 


3 


17.3430879705857 


17.3433 c 



Table III: Energy levels for the x 8 potential from the TBA com- 
pared with previous results: a from table 2 of M], and c from 
table VII of [|, rescaled by 2 4 / 5 

The story might have ended here, but in fact it goes considerably further. Fol- 
lowing ||, we begin by asking about potentials of odd degree, so that the confining 
potential is \x\ 2 , with M now allowed to be a half-integer. It helps to split 
the eigenvalues according to the parity of their eigenfunctions, decomposing D(E) 
accordingly as D(E) = D + (E)D-(E), with 



dhe)=d±(?) n ( i+ f)- 

7 even \ K / 



(16) 



These spectral subdeterminants together satisfy a rather simpler equation than that 
obeyed by the full spectral determinant, which also holds if M is a half-integer Q: 

^^D+iVl-^D-^lE) - VL- 1/2 D + (VLE)D-(Q,- l E) = 2% (17) 

where O = e ln /( M + 1 ). This is very similar to the so-called 'quantum Wronskian' 
condition satisfied by the Q-operators introduced in ref. |21]. The similarity 



becomes more striking when the condition is written in terms of the operators 
A ± (A)^A^ 2 Q ± (A): 

(? 2P ^ 2 A + ( (? 1 /2 A)A _ (g -i/2 A) _ (? -2P/^ A+(g -i/2 A)A _ (? i/2 A) = 2isin(27rP) . (18) 
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The parameters q and /3 2 are related by q = e* 71 "^ 2 . As it stands, this is an operator 
equation, but it becomes a functional equation when applied to the simultane- 
ous eigenvectors \a,p) of the operators A±(A) and P, defined via A±(\)\a,p) = 
A±(X,p)\a,p), P\a,p) = p\a,p). (We refer the reader to |^IJ for the background 
to these definitions.) For brevity we will leave the p-dependence of the eigenvalues 
A±(\,p) implicit from now on; they are entire functions of A 2 , with A±(0) = 1, 
and a finite number of complex and negative real zeroes. The remaining zeroes 
accumulate towards +00 along the positive real axis of the A 2 -plane. To choose 
particular eigenvalues as candidate spectral subdeterminants, we recall that the 
zeroes of D ± (E) lie on the negative real E axis. This selects the 'vacuum eigen- 
values' A~\\), for which all of the zeroes lie on the real axis of the A 2 -plane, and 

suggests to identify A^} '(vE 1 / 2 ) with a ± D ± (— E), using the following dictionary, 
where as before [i = (M+1)/2M : 

f3 2 = 1/(M+1) , p = l/(4M+4) 
;/ = (2M+2)" 1 / 2 'T(^r 1 

a± = v ^(2Af+2)^r(I± 3 L)~ 1 (19) 

Note, a + a~ = sin7r/(2M+2). The constant v is fixed by comparing the behaviour 
of A T (X) as A 2 -» -00 H with that of D ± (E) as E -> +00 |]: 

lnA T (A) ~ ( M +l)r(^) 2 %(-A 2 r ; 

\nD ± (E) ~ \a Q E^ ( a = bo/ (2 sin (m) ) . (20) 

Finally, the zeroes of A^\\) should all lie on the positive real axis of the A 2 -plane 
if they are to map onto those of D^(E). This holds if =F2p > — /? 2 [21], a condition 
which is indeed met here. For a more precise check, we sought some numerical 
evidence. As in plj, consider the functions (A) = e ±imp A^± (qX) / A^± {q^ 1 X) . 
The so-called T-Q relation implies that they assume the value —1 precisely at the 
zeroes either of A± (A), or of a related entire function T(A). For the vacuum 
eigenvalues, the zeroes of T(A) are away from the positive real axis and so a search 
of this line for zeroes of a^/(A)+l will allow us to locate the zeroes of A±\\). At 
the values of p and given by (|l|), the functions f±(0) = Ina^ (e 9 / 2 ^) solve the 
following nonlinear integral equations (NLIE): 



f±(0) = -\ib Q v- 2 ^e e + / p(6-6')ln(l+e f±{e,) )d6' 

JCi 

- [ (f(e-6')ln(l+e- f± ^)d9' ±m/2 (21) 

JCo 
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where the contours C\ and C2 run from —00 to +00, just below and just above the 
real #-axis, and 



e iuje smh% (£-!> du 
2 cosh sinh 2-7T 



(22) 



Such equations first arose in refs. [22, 23|, in the contexts of the (twisted) XXZ 
model, and the sine-Gordon model at coupling /3. 

Solving (pl| ) numerically, we can now test the conjecture (19). For M = 2,3,4, 
the results of Tables I, II, III were reproduced, with disagreements being typically 
in the last quoted digit of the TBA data. Next we set M = 3/2, and obtained 
the results quoted in H for the |x| 3 potential. It was then natural to conjecture 
that the identification remains valid at arbitrary M > 1. In the absence of suitable 
published data, we used the Maple package to diagonalise the Hamiltonian ([!]) in 
a basis of harmonic oscillator eigenfunctions, as in ||. Agreement with ( ^l|) was 
confirmed for various potentials |x| 2M ; some results for M = 15/8 are shown in 
Table IV. 



k 


E k (NLIE) 


E k (QM) 





1.0503451122723 


1.05034511 


1 


3.7190710425856 


3.71907104 


2 


7.2061514537816 


7.20615145 


3 


11.148641889036 


11.1486419 



Table IV: Energy levels for the |a;| 15 / 4 potential computed using 
eq. (f2l|), compared with direct QM results. 



For M < 1, the formulae for the determinants become divergent and need 
further regularisation ||; at the same time, the calculations of |21] depart from 



the so-called 'semiclassical domain' and must be modified. Nevertheless, we have 
evidence that the correspondence continues to hold. At M=l, the sine-Gordon 



model is at the free-fermion point, the kernel (22) vanishes, and the energy levels 
Ek = (2/c+l) of the simple harmonic oscillator are easily recovered. Then at 
M=l/2, the D^(E) are known in closed form Q, leading to the predictions 



» 



(vE 1 ' 2 ) 



nAi(-n 2 E)/Ai(-n- 2 E), 

Q- 1 Ai'(-n 2 E)/Ai'(-n- 2 E) 



(23) 



where Ai(E) is the Airy function and O = e 27r4//3 . These were verified to 15 digits. 
Note that /3 2 =2/3 for M=l/2: this is the A^=2 supersymmetric point of the sine- 
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Gordon model, and it is tempting to conjecture a link with the Painleve III results 



of ref. [24 1, though this remains to be elucidated. Finally, we made a numerical 
check against Maple results at M=7/8, again finding agreement. 

When M is an integer the potential is analytic; it is interesting that these 
cases are mapped onto the reflectionless points of the sine-Gordon model. In the 
TBA framework, these are described by -D2M-i-related systems, with the twist 
p = l/(4M+4) implemented by introducing fugacities ±i on the fork nodes 2M—2 



and 2M—1 (see [24] for similar manoeuvres in the repulsive regime). It can then 
be checked that, for the ground state with Eim-2 = £2M-i, this is equivalent to an 
A2jw_i-related system, thus making a link with the approach described in the first 
half of this paper. 
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